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ABSTRACT 

Using data from the DEEP2 galaxy redshift survey and the All Wavelength Extended Groth Strip Interna¬ 
tional Survey we obtain stacked X-ray maps of galaxies at 0.7 < z < 1.0 as a function of stellar mass. We 
compute the total X-ray counts of these galaxies and show that in the soft band (0.5-2 kev) there exists a signif¬ 
icant correlation between galaxy X-ray counts and stellar mass at these redshifts. The best-fit relation between 
X-ray counts and stellar mass can be characterized by a power law with a slope of 0.58 ± 0.1. We do not find any 
correlation between stellar mass and X-ray luminosities in the hard (2-7 kev) and ultra-hard (4-7kev) bands. 

The derived hardness ratios of our galaxies suggest that the X-ray emission is degenerate between two spectral 
models, namely point-like power-law emission and extended plasma emission in the interstellar medium. This 
is similar to what has been observed in low redshift galaxies. Using a simple spectral model where half of 
the emission comes from power-law sources and the other half from the extended hot halo we derive the X-ray 
luminosities of our galaxies. The soft X-ray luminosities of our galaxies lie in the range 10^®-8 x 10^® ergs s^ . 
Dividing our galaxy sample by the criteria U — B > 1, we find no evidence that our results for X-ray scaling 
relations depend on optical color. 

Subject headings: galaxies; active, galaxies:ISM, ICM, AGNigeneral, XraysiISM, galaxies 


1. INTRODUCTION 

The first systematic studies of X-ray emission in nor¬ 
mal galaxies, conducted using the Einstein Observatory, re¬ 
vealed extend ed and comple x X-ray structure in nearby 
galaxies (e.g.. lEabbianol[T98^ . Subsequent X-ray observa¬ 
tions at higher a ngular resolution with ROSAT augmente d 
these studies (e.g.. lTanaka et alJl994tlMatsushita et alJl994l) . 
but the observations of the Chandra and XMM-Newton 
satellites have revolutionized our knowledge of the X-rav 
source populations withi n galaxies JOwen & Warwick! 120091 : 
iRovilos et al.l l200M and lEabbianol 1201^ for a review). The 
main components contributing to the X-ray emission in 
galaxies can be broadly classified into two types—namely, 
diffuse emission from the hot interstellar medium (ISM) 
which has a cosmological origin, and pointed emission from 
astrophysical objects within galaxies such as X-ray bina- 
ries, evolved stellar populat i ons, and planetary nebula (e.g. , 
Mathews & Brighentil l20()3t iBrassington et alJ l2008l 1200^ 
Bregman & Parriot3l2009HEabbianoll2Q12lk 

The existence of hot gaseous halos in galaxi es is an integral 
part o f the theory of galaxy formation (e.g., IWhite & ReesI 
Il978h . Gaseous halos should be formed when gas condenses 
on to dark matter halos and undergoes subsequent shock heat¬ 
ing and adiabatic compression. Based on this idea, the exis¬ 
tence of exte nded soft X-ray emissi on around galaxies was 
predicted by IWhite & Eren^ dl991h . This established the 


foundation for interpreting galaxy evolution within the Cold 
Dark Matter paradigm. Hot gas halos have been seen in 
early type galaxies since the advent of the Einstein satel¬ 
lite (e.g.. iForman et alJ 1 19851 iTrinchieri & Fabbianol 119851 : 
ICanizares et al.l 1987h and work on the phenomenon has been 


dra ("e.g.-IO’ Sullivan et al. 

l200lUDav 

d et alJl2006l ISun et al. 

2007 

; iJeltema et alJ 12008 

; ISun et al. 

120091 lAnderson et al. 


2013). The detection and characterization of hot halos in 


galaxies has firmly established the basic paradigm that galax¬ 
ies form within collapsing gaseous halos. 

The origin of X-ray emission from astrophysical sources 
within galaxies has mainly been studied thro ugh their cor¬ 
relati o n with galaxy stellar luminos it ies (e.g.. iForman et alJ 
19851 ITrinchieri & Fabbianol 1 1985 1 ICaniz^es et aLl ' 1987 


Sarazin et al.l l2000l iKim & Fabbianol l2003l ISivakoff et alJ 
2003h . Correlations between X-ray properties and stellar lu¬ 


minosity corroborate the connection between galactic X-ray 
emission and the stellar component in galaxies. The rela¬ 
tive fraction of X-ray emission from hot gas in the halo to 
that of the stellar component in the galaxy—in order to sep¬ 
arate astrophysical signals from cosmological sig nals—has 
been addressed in local galaxies with Chandra (e.g.. lFabbianol 
l2006h . We note that, in the context of X-ray star-forming 
galaxies, the diffuse emission is often referred to the X-ray 
emission in the ISM which has been heated by supernovae 
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Fig. 1.— Distribution of stellar mass in each stellar mass bin. The redshifts of these galaxies ai‘e 0.7 < z < 1.0. The redshift distributions are shown in Fig. 2. 
Any galaxy that has been identified as an X-ray bright AGN in the AEGIS-X catalog (L09) has been discarded from our analysis, since we want to minimize the 
X-ray nuclear emission from AGN. 


(e.g., I^ickland et al.ll2002L iRanalli et alj|200^ iMineo et al.l 
l2012h . We clarify that in the current paper, the diffuse emis¬ 
sion refers to the ISM gas that is shock-heated due to its infall 
into the gravitational potential well of the galaxy. 

The results show that the relative fraction of X-ray emission 
from hot halo and stellar components depends on the lu minos- 
ity and morphology of the galaxy ( see lFabbianol200^ and ref¬ 
erences therein). Further correlations appear to exist between 
the X-ray gas ar ound galaxies and galaxy environments on 
large sca les te.g..lWhite & Sarazi^ll99lMBrown & Bregmaij 
200(1 lO’ Sullivan et al.l 1200 It iHelsdon et al.1 UOOlF 


Ellis & O’ Sullivan! l2006t iMulchaev & Jeltemal 1201 


However all of these studies have been carried out for sam¬ 
ples of local galaxies. To check whether the X-ray properties 
of local galaxies are ubiquitous features of galaxy evolution, 
it is important to search for similar correlations at higher 
redshift. 

In this paper we present a study of the relationship between 
X-ray luminosity and stellar mass in galaxies approaching a 
redshift of z ^ 1. Our goal is to characterize the X-ray emis¬ 
sion from galactic sources and compa re them with the re¬ 
sults obtained for loc al galaxies (e.g.. iBoroson et akl 1201 1! 

I Anderson et al.l FoiSh . We perform a stacking analysis of 
galaxies at 0.7 < z < 1.0, observed in the DEEP2 survey 
(iDavis et^l200^ iN^man et al.ll2013h by using X-ray data 
from the All Wa velength Extended Groth Strip Interna tional 
Survey (AEGIS: iDavis et alll2007l: If^dra et al.ll2007h . We 
then compute the X-ray properties of these galaxies from the 
stacked map as a function of galaxy stellar mass and statisti¬ 
cally evaluate our findings in light of local galaxies. AEGIS 
provides a unique combination of wide area and substantial 
depth, and DEEP2 provides a high sampling rate of spec¬ 


troscopy. In combination, these factors allow studies of the 
sources of X-ray emission in galaxies to be extended to higher 
redshift. 

Our paper is organized as follows. In §2 we give a brief de¬ 
scription of our data sets and describe our methodology. We 
present and discuss our results in §3 and §4 respectively. We 
finally summarize our conclusions in §5. Throughout the pa¬ 
per we assume a spatially flat, ACDM cosmology: = 0 .28, 

Oa = 0.72, Q.b = 0.04, and h= 0.71 (ISpergel et al.ll200% . 

2. DATA SETS 

We use data from the AEGIS-X survey (iLaird et akl120091 : 
L09 hereafter) and the DEEP2 galaxy redshift survey. In this 
section we will provide a brief description of our data sets. 

The DEEP2 Redshift S urvey (Deep Evolutionary Ex¬ 
ploratory Probe; IDavis et ani2003h is a project that has mea¬ 
sured redshi fts of 50,()00 gala xies to a lir mting magnitude of 
Rab = 24.1 dCoil et al.ll20^|Paber et 101200% . The DEEP2 
Galaxy Redshift Survey is a spectr oscopic survey of galaxies 
using the DEIMOS spectrograph dPaber et alJl2003h on the 
10 m Keck II telescope to survey optical galaxies at z ^ 1 in 
a comoving volume of approximately 5 x lO^h^^Mpc^. The 
survey covers 3 deg^, over four widely separated fields to mit¬ 
igate cosmic variance. Due to the high resolution (R ^ 5000) 
of the DEEP2 spectra, redshift errors, are ~ 30kms^^ For 
more details of DEEP 2 observations, catalog c o nstruction, 
and data red uction see IDavis et alJ d2003L l2005h : ICoil et af] 
d2004bh . and iFaber et al.l d2007h . 

To measure accurate stellar masses, it is necessary to 
have photometry extending from the visible into the near- 
infrai'ed (4,000-22,000 A). However, since only ^ 33% 
of DEEP2 galaxies with redshifts have stellar mass mea- 
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Fig. 2.— Distribution of redshift in each stellar mass bin. §2 details the construction of samples with matched redshift distributions. The stellar masses are 
shown in each panel. 


surements from K band magnitudes (iBundv et al.l 120061: 
IConselice et alJl^OSh . we use a different method to measure 
the stellar masses of the remaining galaxies. _ 

We used the method proposed bv iLin et^ (l2007h and more 
recently bv I Weiner et al.l (l2009h to estimate stellar masses for 
the remainder of our galaxy sample. The rest frame Mb and 
U — B color* are us ed to estimate the ste l lar masses using 
color-M/L relations (iBell & de Jongll^OOlt iBell et af] 120031) 
an d then corrected usin g X-band photometry and SED fitting 
bv iBundv et al.l (l2006ll . We have 14121 galaxies in the Ex¬ 
tended Groth Strip (EGS) field that have measurements of 
stellar mass. We apply a redshift cut (0.7 < z < 1.0) to our 
parent galaxy sample. We also discard galaxies from our sam¬ 
ple that are identified as X-ray bright AGN in L09. Our final 
galaxy sample comprises 2914 galaxies. We subdivide our 
galaxy sample into six stellar mass bins shown in Eig. 1. 

In Eig. 2 we present the redshift distribution of our sources. 
We expect that galaxy X-ray properties will depend on 
redshift—since galaxy ages, luminosity distance, angular di¬ 
ameter distance, and surface brightness dimming will all 
evolve with redshift. Hence, to compare the luminosity de¬ 
pendence with galaxy stellar mass we require the redshift dis¬ 
tributions to be identical in each stellar mass bin (shown in 
Eig. 1). We note that the redshift distribution of the two low¬ 
est stellar mass bins (top-left and top-middle panels in Eig. 2; 
the figure actually shows the matched distribution described 
below) differ from the other four mass bins. 

To ensure that we can study low and high stellar mass sam¬ 
ples over a similar volume, we applied a redshift matching 
criterion. To obtain matched redshift distributions we stack 

* We obtain rest frame (U-B) colors using the method of IWillmer et alj 

f200a . 


some galaxies multiple times for the lower stellar mass bins. 
We note that a unique sample of galaxies do not reproduce the 
redshift distribution in the two lowest stellar mass bins. Thus, 
it is necessary to repeat objects in order to match the distribu¬ 
tion of our samples. The final redshift distributions for each 
stellar mass bin are shown in Eig. 2. We note that this tech¬ 
nique could be problematic in cases where a single object may 
come to dominate the properties of the stellar mass bin (e.g., 
if it is the brightest). 

The AEGIS-X survey consists of 8 deep Chandra ACIS- 
I fields, each with a total integration time of about 200ks. 
The total area covered in the survey is about 0.67deg^. The 
images are constructed in four energy bands namely, 0.5- 
7.0keV (full), 0.5-2.0keV (soft), 2.0-7.0keV (hard) and 4.0- 
7.0keV (ultra-hard). The limiting fluxes are 2.37 x 10^*®, 
5.31 X 10-*^ 3.76 X 10-*^ and 6.24 x lO^'^ergs-* cm-2, 
respectively. The limiting flux is defined as the flux to which 
at least 1% of the survey area is sensitive. ^ A point source 
catalog in the EGS held has been provided in L09. The cat¬ 
alog consists of a total of 1325 sources. We identified X-ray 
bright point sources in our galaxy catalog from this parent cat¬ 
alog. We also used the AEGIS-X point source catalog to mask 
out all X-ray bright point sources in the EGS held. Note that 
we also discarded X-ray bright AGN from our galaxy sample. 
The X-ray bright AGN were identified from the L09 point 
source catalog. We used a luminosity cut (Lx > lO'^'ergs^*) 
to identify the X-ray bright AGN from our source catalog. 

2.1. Methodology 

We perform a stacking analysis to identify the X-ray emis¬ 
sion, from our galaxy sources. The event files, the expo- 

^ See L09 for more extensive details of how the data were reduced 
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Fig. 3.— Correlation between stellar mass and X-ray counts in the three energy bands. The left-hand, middle, and right-hand panels represent the soft, hard 
and ultra-hard bands. The data points with error bars are detections whereas the open circles (with arrows) are 2o upper limits. The dashed line in each panel 
represents the best-fit power law. The best-fit slopes (a = dlogAl/dlogM*: N and M* being X-ray luminosity and stellar mass respectively) are shown in the each 
panel. When converted to X-ray luminosity (Fig. 5), the best-fit slope a = dlogL/dlogM* is exactly identical to what has been observed with photon counts. We 
have included both the significant and the insignificant detections for obtaining our best-fits. We note that the slopes for the soft and the hard bands are similar 
while the best-fit slope for the ultra-hard band is consistent with zero. We also note that the correlation co-efficient is only significant for the soft band. The errors 
on the stellar masses are the 1 a spread in the distributions and are identical in each panel. See §4 for more details. 


sure time maps and the effective area maps are provided in 
L09. We first correct our event files using the exposure time- 
effective area maps. We call these effective area-exposure 
corrected (EAEC hereafter) maps. The gaps and singularities 
in the EAEC maps are corrected by assigning zero counts to 
those pixels where the EAEC maps have singularities (due to 
zero exposure and/or effective area). We then identify sources 
in the EAEC maps and select a 5 x 5 square arcminute region 
around each source and sum them to construct the stacked im¬ 
age. 

We identify the point sources in each map from the point 
source catalog presented in L09. We calculate the point 
spread function (PSE) for each point source using the tech¬ 
nique described in L09. The gefps/routine, provided by L09, 
is used to obtain the PSEs. We then mask the point sources us¬ 
ing the ellipse corresponding to the 95% encircled energy ra¬ 
dius (EER). We mask the region which encompasses a circu¬ 
lar area with a radius 1.5 times that of the semi-major axis of 
the 95% encircled energy ellipse of a particular point source. 
This allows us to create a conservative mask for each point 
source. 

We also construct a stacked mask map. It is important to 
note that a source in one map could be a potential contaminant 
point source (to be masked) in another map. The EAEC maps 
are divided by the stacked mask maps to obtain the average 
photon counts in each pixel of our final stacked maps. We 
construct a total of 18 maps for the three bands. In each band 
we have six stacked maps corresponding to the six stellar mass 
bins shown in Eig. 1. 

To calculate the background, we extract the counts from the 
stacked maps in an annular region within 50-70" of the cen¬ 
tral source. Since the PSE of the central sources are ^ 2-3" 
the region for background extraction is, by construction, at 
a substantial distance from the sources. We find that chang¬ 
ing the area of the region for background extraction does not 
affect the measured background as long as the region is suf¬ 
ficiently far from the sources. The background counts are 
estimated to be 1.68 x 10^'*’ cm^^s^'pixel^^ 3.4 x 10^'** 
cm^^s^'pixel^^ and 2.2 x lO^**' cm^^s^^pixeR* respec¬ 
tively for the soft, hard, and ultra-hard bands. 

3. RESULTS 

Eig. 3 shows the total X-ray counts of our galaxies as a 
function of stellar mass. The left-hand, middle and right-hand 
panels are for the soft, hard, and ultra-hard bands. The total 


counts are calculated using an extraction scale of 3 arcsec- 
onds. Note that the redshifts of these sources are similar and 
hence the total counts are fair representations of their lumi¬ 
nosities. To estimate errors we compute the standard devia¬ 
tion of counts within the extraction region and then divide by 
the square-root of the total number of pixels in the circle to 
obtain the standard error on the mean. Since we are interested 
in the total counts, we multiply the standard error by the to¬ 
tal number of pixels. Hence the error-bars are O'/N, where o 
and N refer to the standard deviation of counts and the total 
number of pixels in the circle. 

It is important to note that according to this definition of 
our error on the counts, the signal-to-noise is proportional to 
\/N. Hence a larger number of pixels within the extraction 
radius would result in a more accurate detection, if all other 
quantities stay the same. The open circles with arrows repre¬ 
sent 2a upper limits of non-detected emission (i.e., the signal 
for these data points is consistent with the background to 2a). 
The dashed line in each panel refers to the best-fit power law. 
We h ave used the weighted l east square fitting method in this 
case (lAndreon & Hurnll2012h . The missing data points depict 
non-detections where the 2o upper limits are either below the 
range shown in Eig. 3 or are consistent with zero counts. 

We also note that in the current method of background ex¬ 
traction the actual amount of background remains uncertain 
since we do not know if a single photon came from the source 
or the background. So an uncertainty on the background 
should be considered while computing the error on the total 
counts from the source. As mentioned above we have ex¬ 
tracted our background from a much larger region compared 
to the source extraction radius. In addition, we have checked 
that the background remains the same if we increase the size 
of the annulus. Due to the large area of the annulus from 
which we have extracted the background we note that the un¬ 
certainty on the background counts will be smaller than the 
Poisson noise we have considered. Hence we did not include 
the associated error in our error budget. 

In the soft band we see a clear correlation between X-ray 
counts and stellar mass. In the hard band (panel b of Eig. 3) 
we do not detect significant X-ray counts for the lower stellar 
mass bins. In the ultra-hard band (panel c of Eig. 3) our de¬ 
tections are less significant than in other bands. We used both 
significant and insignificant detections while computing the 
best-fit power-laws for Eig. 3. We note that the soft and the 
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hard bands reflect similar slopes for the power-law. The slope 
in the ultra-hard band is consistent with zero. It is, thus likely 
that the bulk of the emission in the soft and the hard bands are 
coming from similar source population. While computing the 
total counts in Fig. 3 we have fixed the extraction radius for 
each stellar mass bin. However we note that since halo masses 
are correlated with stellar masses the sizes of galaxies and the 
corresponding dark matter halos will be different in each case. 
We thus compute the counts using scales corresponding to a 
fixed fraction (/) of R 200 (defined as the radius at which the 
mean density of the halo is 200 times the critical density of 
the Universe) derived from the respective halo masses. R 200 
is computed from the dark matter masse s using the theoretica l 
scaling relation with stellar mass from iMoster et al.l (l2010h . 
We note that the best-fit slopes are statistically identical in 
each case. 

From the counts in Fig. 3 we compute the hardness ratios. 
The hardness ratio is defined as {H — S)/{H + S), where H 
corresponds to the total counts in the hard band and S refers 
to the total counts in the soft band. The hardness ratios are 
shown in Fig. 4. Note that for the non-detections we use 
the observed value of the counts and their observed errors. 
This accounts for the large error-bars in the hardness ratios in 
the lower stellar mass bins. To compare our hardness ratios 
with different spectral models of emission we compute model 
hardness ratios for a fixed flux (10^^® ergs s^* at soft X-ray 
energies) with different spectral models using the web-based 
PIMMS tool.3 

The different models are a) power law with positive pho¬ 
ton index, b) thermal black body, c) Bremsstrah lung, and d) 
the A strophysical Plasma Emission Code (APEC: ISmith et all 
I2001h. We assuine a G alactic Nh of 1.3 x lO^'^cm ^ 
(iDickev & LockmunI Il990l) and a redshift of z = 0.8. To 
study the effect of intrinsic absorption on the hardness ra¬ 
tios we calculated the model hardness ratios using an intrinsic 
Nh — 5 X lO^^^’ergs^^ We note that the change in hardness 
ratios due to absorption in some models is insignificant for the 
purpose of our analysis. 

We And that in our definition of the hardness ratio, the 
pure black-body spectrum gives positive hardness ratios for 
KbT > l.Skev. Our data thus exclude the pure black- 
body model for these energies. We compute the hardness 
ratios from the APEC model using different metallicities 
since the plasma e mission will depend on metal abundances 
(ISmith et alJl2001h . We And that the differences in the hard¬ 
ness ratios due to differences in metallicities are well within 
our error-bars. So for our purpose we compute hardness ra¬ 
tios using solar abundances. In the top panel of Eig. 4 we plot 
the hardness ratios for the APEC model with solar metallicity 
as a function of plasma temperature. The Brehmsstrahlung 
model is depicted by blue triangles, as a function of tempera¬ 
ture. In the bottom panel we plot the derived hardness ratios 
for a power-law model, as a function of photon index. The 
black squares correspond to the measured hardness ratios of 
our galaxy sources. 

We note that X-ray emission from different sources will 
originate from different spatial scales in galaxy X-ray profiles, 
and ideally both spectral and spatial resolution of the sources 
will give us information about the exact source population. 
However, we do not have the spectral or spatial capability 
to distinguish the point sources in our high redshift galax- 
ies from the extended emission. Recently, lAnderson et alJ 

^ http://heasarc.nasa.gov/Tools/w3pimms.html 


(1201 3h performed a detailed stacking analysis of nearby iso¬ 
lated galaxies using ROSAT All Sky Survey data (IVoges et alJ 
Il999ll to search for hot halos. At redshifts below 0.06 their re¬ 
sults show that the X-ray emission in isolated galaxies within 
a spatial scale of 50 kpc from the center of the galaxy is very 
likely to be a linear combination of a power-law emission 
from X-ray binaries and an APEC plasma emission from the 
h ot halo. _ 

lAnderson et alJ (1201 3h performed detailed simulations to 
identify the degree of this degeneracy in isolated local galax¬ 
ies and showed that the exact fraction of extended emission 
to power-law emission depends on the spatial scale and other 
parameters of the galaxy. We thus follow a simpler approach 
to compute the luminosities of our source galaxies. 

Eollowing the arguments of lAnderson et alJ (1201 3h we ex¬ 
tract the X-ray luminosities of our galaxies assuming an equal 
contribution from both extended emission and power-law 
emission. Our measured hard ness ratios also suggest a similar 
level of degeneracy as seen bv lAnderson et alJ (120131) . mostly 
due to the errors in our measured hardness ratios.^ We note 
that the assumption of Gaussianity in the error-propagation 
for small number counts can lead to a 30% overestimate of 
our errors. A more accurate method involves the Bayesian 
Estimation of Hardness Ratios (iPark et alJl2006h approach. 

4. DISCUSSIONS 

Previous studies have investigated the contribution from 
different kinds of sources to the X-ray luminosity of nearby 
galaxies, and their respective correlation with the stellar mass 
or the X-band o r the B-band luminosity o f the galaxy (e.g.. 


GilfanovI 120041 iBogdan & GilfanovI 120081 iRevnivtsev et al. 

2008HBogdan & Gilfanov 

201 lUBoroson et alJ201 lb. Recent 

work (iBoroson et al.ll201 1 
correlation between X-ra} 
in a galaxy varies for difl 
ferent point sources are ste 
unresolved X-ray binaries 
ray binaries (LMXB). Th 
high redshift galaxies will 
of these sources. Also 01 
contribution from unresol 
dividually detected in our 
to the signal in the co-adc 
the typical photon index ft 
2.5. This spectral index 

, B11 hereafter) has shown that the 
/ luminosity and stellar population 
erent source populations. The dif- 
fllar X-ray sources (ABs and CVs), 
(XRB), and resolved low mass X- 
us the power-law emission in our 
be the integrated emission from all 
IT data will possibly include some 
ved AGN, which have not been in- 
X-ray catalogs but can contribute 
ed maps. Prom Pig. 4 we And that 
31' our sources lies in the range 1.5- 
falls in the range of both LMXBs 

(e. 2 .. iGrimm et al.ll2002l; lPabbianoll2006l; 

Burke et al.ll2013b. 

and AGN sources (e.g., |f 

201 iHLanzuisi et alJl2013 

vlandra & Poundslll994l; iScott et al. 
). 

AGN activity is mostly seen in 
dback from AGN has been linked 
formation fe.g.. iKauffmann et al. 

It is well known that 
massive galaxies and fee 
with suppression of star 

20031 fDunloD et al.l 120031: ISchawinski et al. 2007: 

Alonso-Herrero et alJ l2008l iBrusa et alJ l2009l IVitale et al. 


2OI3I) . but the exact correlation between AGN feedback and 


galaxy stellar mass is yet to be confirmed. Some studies 
suggest that there i s no correlation b etween stellar mass and 
AGN activity (e.g.. lAird et al.liy012l) . Hence we assume that 
although nuclear X-ray emission from unresolved AGN can 
contaminate the X-ray emission observed at the high stellar 
mass end of Pig. 3 (see panel a), it will not be responsible 

^ Note that our deri ved model hardness ratios are slightly different from 
lAnderson et alj 120131) due to differences in the telescope response function 
and exact energy scales of the soft and hard bands. 
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Fig. 4.— Hardness ratios for different spectral models. The top panel 
shows the hardness ratios for the Bremsstrahlung model (blue solid line with 
triangles), and the APEC plasma model as a function of plasma tempera¬ 
ture. In the bottom panel the red solid line with open circles shows the 
value of hardness ratios for a power law model with different photon indices. 
The black squares in each panel represent the measured hardness ratios of 
our sample as a function of stellar mass. The hardness ratios for the non¬ 
detections are calculated using the exact observed value of the counts and the 
exact observed error. That explains the large eiror-bars in the lower stellar 
mass bins. The hardness ratios seem to be degenerate between power law 
emission and APEC/Brem sstrahlung emission as is also observed in the case 
of low redshift galaxies bv lAnderson et alj 120131) . See §4 for more details. 


for the emission observed in the entire stellar mass range. 
Thus we adopt a photon index of 1.5, typically assumed for 
emission from stellar sources in galaxies. We note that the 
supernova heated X-ray plasma in star-forming galaxies can 
have photon indices of 2. We stress that our data does not 
exclude such possibilities. However for the purpose of our 
luminosity calculations we note that a photon index of 1.5 
and 2.0 does not affect our results signihcantly. We thus 
stress that our results could also represent the plasma in 
star-forming galaxies. 

To understand the ISM emission in our galaxies we need to 
have knowledge abo ut their host dark matter halos. Recently 
iMostek et al.l ( 20131) studied the clustering of DEEP2 galaxies 
at z ~ 1 as a function of stellar mass, hnding that clustering 
strongly depends on stellar mass for blue star-forming galax¬ 
ies in the mass range 9.5 < log(M*/M 0 ) <11, but has no 
dependence on stellar mass for red quiescent galaxies. The 
host dark matter halos of the blue galaxies are in the range 
12.32 < logMh,io[h-'M 0 ] < 12.75, whereas red galaxies in 
the stellar mass range 10.5 < log(M*/M 0 ) < 11.5 are re¬ 
siding in host halos of masses 12.95 < logMjjjj[Q[h^^M 0 ] < 
13.05. In summary, clustering results suggest that galaxies 



Fig. 5.— Lx —Lb relation for our high redshift galaxies. The X-ray lu¬ 
minosity has been calculated using a simple two-component spectral model 
where 50% of the emission comes from the stellar sources in the galaxy 
(power-law emission) and the other 50% is assumed to come from the hot 
ISM (APEC plasma emission; see Fig. 4 for the degeneracies in the hardness 
ratios between different spectral models). The spread in fi-band luminos¬ 
ity is large due to our primai'ily mass-selected sample. The best-fit power- 
law slope is 1.25 ± 0.5, consi stent with the results obtained for local galax- 
ies using ^-band lum inosities IHelsdon et alJ2QCT) . and X-band luminosities 
iBoroson et alj20ll1) . 


with stellar mass > lO**'h^^M 0 generally reside in massive 
galaxy/galaxy group scale halos. 

To check the distribution of galaxy colors in our sample we 
apply a high redshift color-cut, where we identify red galax¬ 
ies with restframe color {U — B) > 1. We find that the frac¬ 
tion of red galaxies varies strongly with stellar mass. This, 
of course, is reflective of the selection criteria for DEEP2, 
where low mass galaxies are selected if they are bright enough 
in B-band. We And that roughly 30-40% of our sources in 
the highest stellar mass bins are red galaxies. Thus, while 
computing the plasma emission for our galaxies, we adopt 
temperatures corresponding to massive galaxy scale halos 
1 O*^'^M 0 ) for our low mass galaxies and group scale halos 
1 O^^M 0 ) for high mass galaxies. Prom the two-component 
fitting (power-law emission with a photon index 1.5 plus a 
massive galaxy/galaxy group scale halo plasma emission) we 
obtain the soft X-ray fluxes of our galaxies and use the me¬ 
dian redshift to convert to soft X-ray luminosities. The results 
are shown in Pig. 5. 

We note that we are interested in the correlation between 
galaxy stellar mass and X-ray luminosity, but the majority of 
the low redshift studies use the X-band or the B-band lumi¬ 
nosity of the galaxy to search for such correlations. We thus 
compute the average B-band luminosities of our stellar mass- 
selected sample and compare our results with local galaxies. 
We emphasize that the X-band luminosity is better correlated 
with stellar mass than B-band luminosity but our restframe X- 
band lumino sities are incomple te mostly due to confusion in 
Spitzer data (iBundv et al.ll2006h . So, we rely on B-band lumi¬ 
nosities. We note that our galaxies are selected based on stel¬ 
lar mass and thus in each mass bin we do observe a substantial 
scatter in B band luminosity. This should be taken into con¬ 
sideration while assessing the strength of the reported corre¬ 
lation. Another important effect that can alter the correlation 
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is the morphological types of our galaxies (most of the low 
redshift scaling relations hold true for sample of early-type 
galaxies). We do not have information on galaxy morphology 
and hence we performed a spectral cut and re-examined our 
scaling relations based on this color-cut. As noted previously, 
the fraction of red galaxies varies strongly with stellar mass 
due to the selection criteria of the DEEP2 survey—neglecting 
that selection effect, the X-ray scaling relations for red galax¬ 
ies are statistically identical to the X-ray scaling relations for 
the full sample. 

We now compare our results with the scaling relations for 
local galaxies. Bll showed that the correlation between X-ray 
luminosity and stellar X-band luminosity in a galaxy varies 
for different source populations. We do not have the spectral 
or spatial resolution to individually resolve the components 
that are responsible for the source population. We will thus 
compare our total X-ray luminosity (which we have obtained 
assuming an equal contribution from the power-law emission 
and extended emission) slopes with Bll and other studies that 
used local samples. Bll argue that for LMXBs the slope 
of the power law is ^ 1 as found in previous studies (e.g . 

IWhite et al.ll2()()l iKim & Eabbi^lTOOl iDavid et al.l[2(m . 

Eor the diffuse hot gas the slope is steeper compared to the 
LMXBs (Bll). There has been debate in the literature re¬ 
garding the scatter of the Lx — Lk relation for hot gas. This 
difference in slope could arise due to differences in the dark 
matter potential well (or the environment) of t he galaxy (e.g. 
iJeltema et al.l2008l;lMulchaev & Jeltemal2010b . Bll attribute 
another component as a contributor to the X-ray emission 
in galaxies; namely unresolved LMXBs (designated “unre¬ 
solved nuclear emission” in Bll). Eor this particular compo¬ 
nent B11 do not find any correlation between X-ray luminos¬ 
ity and X-band luminosity. 

The Lx — Lk correlations for different X-ray sources in low 
redshift galaxies are summarized in Eig. 2 of B11. Similar re¬ 
sults should also be expected for B-band luminosities (Bll). 
Our best-fit slope for the Lx — Lb scaling shown in Eig. 5 is 
1.25 ±0.5, consistent with the total Lx — Lk slope o f B11 
(1.4 ± 0.2) and the Lx — Lb slope of iHelsdon et al.l (120011: 
1.5 ± 0.2) and lShaplev et al.l (1200 it 1.5 ±0.1). Although this 
is consistent with our findings we emphasize that our scaling 
relations with luminosity are likely to be affected by both se¬ 
lection biases and theoretical biases (e.g. assumption of the 
nature of the source populations) as discussed previously. 

5. SUMMARY 


In this study we have detected significant X-ray emission 
from normal galaxies at high redshift (^0.8) as a function of 
stellar mass. We have shown that there exists a significant 
correlation between stellar mass and X-ray emission at scales 
within lOkpc in the soft X-ray band in these high redshift 
galaxies. We derived a power-law relationship between in¬ 
tegrated X-ray counts and galaxy stellar mass. The slope of 
the best-fit power law is 0.58 ±0.1. We compute the hard¬ 
ness ratios of our sources and show that the hardness ratios 
of our sources are degenerate between two spectral models, 
namely point-like power-law emission and extended plasma 
emission in the ISM. T his has been robustly demonstrated by 
lAnderson et al.l (1201 3h in the case of nearby galaxies and our 
results are similar to those observations. 

Using a simple two-component spectral model where half 
of the emission comes from hot halos and the other half comes 
from stellar point sources in the galaxy, we compute the ap¬ 
proximate luminosities (10^® < < 10^* ergs s^') of our 

sources. We also compute the B-band luminosities of our 
sources and show that the Lx — Lb scaling relation is con¬ 
sistent with low-redshift results (Bll). We also made a color 
cut (t/ — B) > 1 to our galaxy sample to examine the effect 
of galaxy color (broadly equivalent with morphology) on the 
scaling relations. The scaling relations are statistically identi¬ 
cal when such a color-cut is imposed, within the limitations of 
the DEEP2 survey selection. Deeper extragalactic X-ray sur¬ 
veys should provide additional insights into the nature of the 
X-ray source populations in high redshift galaxies, which can 
reveal important clues about galaxy evolution. Our analysis 
provides a first step in this direction. 
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